Directed translocation of a flexible polymer through a cone-shaped channel 
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Entropy-driven directed translocation of a flexible polymer through a cone-shaped channel is 
studied theoretically and using computer simulation. For a given length of the channel, the effective 
force of entropic origin acting on the polymer is calculated as a function of the apex angle of the 
channel. It is found that the translocation time is a non-monotonic function of the apex angle. By 
increasing the apex angle from zero, the translocation time shows a minimum and then a maximum. 
Also, it is found that regardless of the value of the apex angle, the translocation time is a uniformly 
decreasing function of the channel length. The results of the theory and the simulation are in good 
qualitative agreement. 

PACS numbers: 87.15. A-, 36.20.Ey, 87.15.H- 



I. INTRODUCTION 



known that the Brownian motion of particles 
presence of asymmetric structures and in non- 
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equilibrium conditions may result in directed motion. 
This phenomenon, called Brownian motor or Brownian 
ratchet is ubiquitous in the living cells. The known ex- 
amples of this phenomenon are the electric potential dif- 
ference through the ion channels and the movement of 
the kinesin motor protein along the microtubule 0, Hj]. 
This phenomenon has attracted great interest in recent 
years, due to its applications in the separation of parti- 
cles 0, HI and making pumps [|| and motors in fine 
dimensions. The directed motion of particles arisen by an 
asymmetric structure as a non-equilibrium phenomenon 
has been observed in systems covering a broad range of 
scales such as macroscopic elastic discs J8|, mesoscopic 
gears microscopic colloidal systems [j], moving cells 
3] and ions @. 

In the case of polymers, directed translocation of a 
polymer through a curved bilayer membrane Q and poly- 
mer passage through a membrane in the presence of chap- 
erons [3, [T3| have been studied. Despite extended stud- 
ies in this field, the polymer motion through asymmet- 
ric structures such as a cone-shaped channel is poorly 
understood. It has been shown that cone-shaped chan- 
nels have important applications in rectifying the ionic 
currents and in the simulation of biological ion channels 
Q . Polymer translocation through cone-shaped channels 



has also been studied, in the literature [11[. It has been 
observed that the dependencies of the translocation time 
and the capture rate on the applied voltage and the poly- 
mer length are qualitatively similar to those of the cylin- 
drical channels. The most important point in the case 
of cone-shaped channels is the very high strength of the 
electric field in the narrow entry of the channel. It causes 
the ionic current through the channel to be affected only 
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by the few monomers in the narrow entry of the channel 
j 1 11 ] . This point is important in the application of these 
channels in DNA sequencing. One important challenge 
of DNA sequencing by the ordinary channels such as a- 
hemolysin protein channel is the simultaneous effect of 
10-15 nucleic acids inside the channel on the ionic cur- 
rent [121 ] . However, in the cone-shaped protein channel, 
MspA, only two nucleic acids which are close to the chan- 
nel apex simultaneously affect the current [13j. In DNA 
translocation experiments through MspA, the difference 
between the effects of the four types of the nucleic acids 
on the ionic current is larger, compared to the experi- 
ments using the cylindrical channels [l3[. This is also 
an important advantage of the cone-shaped channels for 
DNA sequencing (l4| . 

In this paper, the translocation of a flexible polymer 
through a cone-shaped channel in the presence of no ex- 
ternal driving field is studied, theoretically and by com- 
puter simulation. During the translocation process, a 
force of entropic origin acts on the part of the polymer 
which is inside the channel. This force originates from the 
entropic tendency of the polymer toward the larger entry 
of the channel. The translocation time is a decreasing 
function of this force. We set out to obtain the effec- 
tive force, /, acting on the polymer as a function of the 
channel apex angle and the channel length and compare 
it with the simulation results. For a given length of the 
cone-shaped channel, we calculate the effective force for 
the two cases of small and large apex angles of the chan- 
nel, theoretically. During the translocation of a polymer 
through a narrow channel in a wall (from the cis to the 
trans side) , it is known that the monomers of the poly- 
mer segment passed to the trans side accumulate near the 
wall, in front of the channel exit. Therefore, to calculate 
the effective force acting on a polymer which translocates 
through a cone-shaped channel of small base diameter 
(from the apex to the base), we consider the channel as a 
closed cavity confining a segment of the polymer. In the 
case of large base diameters, however, we consider one 
of the polymer ends fixed in the apex of the cone-shaped 
channel and obtain the confinement free-energy and then 
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the exerted force to the polymer. As a further check of 
the reliability of the results, the force is calculated from 
two different methods in the case of large diameters of 
the channel base. 

For a given length of the channel, we find that the ef- 
fective force is a non-monotonic function of the channel 
apex angle. Combination of the results obtained from the 
two cases of small and large apex angles shows that the 
force as a function of the apex angle has a maximum and 
then a minimum. We also obtain the dependence of the 
force on the channel length. For each value of the chan- 
nel apex angle, the force is a monotonically increasing 
function of the channel length. These results are sup- 
ported by our simulation data for the translocation time 
of a polymer through a cone-shaped channel. The the- 
ory also predicts the polymer length inside the channel, 
which is in agreement with the simulation results. 

The rest of the paper is organized as follows. In Sec. UH 
we present the theory, for the two cases of small and large 
apex angles. The force on the polymer and the average 
number of the monomers inside the channel are discussed. 
In Sec. IIII1 the simulation results of the polymer translo- 
cation through a cone-shaped channel are presented and 
compared with the theory. Finally, in the last section, the 
paper is summarized, the results are discussed and some 
notes are presented on the polymer equilibrium during 
its translocation through the cone-shaped channel. 



II. THEORY 

Simulation results (in the next section) show that the 
polymer translocation through a cone-shaped channel is 
a driven process. When one end of a polymer is fixed 
in the channel apex, a driving force is exerted on the 
polymer, by the channel. This force results from the in- 
crease in the polymer entropy with moving to the wider 
parts of the channel. During the polymer translocation, 
also, the very narrow apex of the channel divides the 
polymer into separate parts with separate free energies. 
The free energy of the polymer part inside the channel 
decreases with the polymer movement toward the wide 
entry of the channel. Indeed, the system is in nonequi- 
librium and writing a total free energy for the polymer is 
meaningless. Hence, in this section, we consider the part 
of the polymer which is inside the channel as a separate 
polymer in a conical channel and calculate its free energy. 
The derivative of the free energy gives the force exerted 
on the polymer, which is the driving force in the polymer 
translocation. We discuss the assumption of polymer in 
equilibrium in Sec. IIVI 

It has already been shown that during the polymer 
translocation through a narrow channel, the monomers 
of the polymer segment passed to the opposite side crowd 
close to the channel [ill [ll| . Accordingly, for the case of 
small apex angles of the channel, we assume that the 
channel is a closed volume containing a part of the poly- 
mer. With this assumption, we calculate the force act- 
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FIG. 1: (Color online) (a) The effective force exerted to a 
polymer confined in a closed frustum of length L = 156, as a 
function of the cone apex angle. As can be seen, the effective 
force has a maximum at small apex angles. Inset: Schematic 
of a polymer confined in a frustum by two walls perpendicular 
to the channel axis, (b) The effective force as a function of 
the frustum length, for a channel of apex angle a — 3 ° . For a 
given value of the apex angle, the effective force monotonically 
increases with the frustum length. Do = 2& for both panels 
(a) and (b). 



ing on the polymer due to the asymmetric shape of the 
channel. In the case of large apex angles, we calculate 
the force acting on a polymer whose end is fixed in the 
apex of a cone-shaped channel and calculate the driving 
force on the polymer, using two different methods. By 
combining the results of the two cases, we conclude the 
qualitative behavior of the force acting on the polymer 
as a function of the channel apex angle. 



A. The case of small apex angles 

Consider a flexible polymer, consisting of N spheri- 
cal monomers of diameter 6, which is confined inside 
a frustum (see the inset of Fig. Ula)). The confine- 
ment free energy of a polymer inside a closed space 

of volume il is known to be F ~ ksT Q b ^ 3 " 1 , 

where v is the Flory exponent fl7^ . In our 
case, the volume of the confining space is equal to 
D, ~ ((D + 2(a + L)tana) 3 - (D + 2a tan a) 3 ). 
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Therefore, the free energy of confinement can be writ- ten as 



F~k R T 



N 3u b 3 tana 



(D + 2(o + L) tan a) 3 - (D + 2a tan a) 2 



1/(3^-1) 



(1) 



r 



Here, a is half of the apex angle of the channel, and a 
is the distance between the region in which the polymer 
is confined and a cross section of the cone with diameter 
Da- L is the length of the region that confines the poly- 
mer. When the confinement region is moved toward the 
wider parts of the channel (increasing a in the inset of 
Fig. HI a)), the entropy of the polymer increases and its 
confinement free-energy decreases. In other words, if one 
removes the two confining walls, which are perpendicular 
to the channel axis, the polymer moves to the wider part 
of the channel to gain more entropy. This is the origin 



of the force exerted to the segment of a polymer inside 
the channel, in the course of its translocation through a 
cone-shaped channel. To obtain the driving force in the 
polymer translocation through a cone-shaped channel of 
length L and tip diameter Do, we calculate the derivative 
of the confinement free-energy (Eq. [1]) with respect to a, 
at a = 0; / = — (^r) a _ - O ne should note here that 
for a translocating polymer, N in Eq. [T]is the number of 
monomers inside the channel, not the total length of the 
polymer. The force, /, acting on the polymer segment 
inside the channel is obtained as 



Do 



tana 



3v 



1 + 2-^ tana) -1 



1 + 2-^- tana ) I 
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In this equation, the number of monomers inside the 
channel, N, should be substituted from Eq. 0] (see be- 
low). In Figs. QJa) andQJb), the force, /, is shown as a 
function of the channel apex angle, a, for a given length 
of the channel, and as a function of the channel length, 
L, for a given value of a. As can be seen in Fig. HJa), 
the force as a function of the apex angle has a maximum 
at small values of a. The existence of this maximum 
can be explained regarding the two determining factors 
in the force: the strength of the polymer confinement in 
the channel and the magnitude of the asymmetry in the 
channel shape. At small apex angles, the channel volume 
is small and the strength of the confinement is high. By 
increasing the apex angle and hence the asymmetry of 
the channel shape, the force of entropic origin exerted 
to the polymer increases. By more increasing the apex 
angle, the confinement effect of the channel on the poly- 
mer weakens and the strength of the force decreases. The 
maximum value of the force corresponds to a value of the 



apex angle for which the combination of the confinement 
effect and the asymmetric shape of the channel has the 
optimum driving effect. 

B. The case of large apex angles 

In the case of large apex angles, considering the poly- 
mer as a confined one in a closed volume is not reason- 
able. Instead, we consider a polymer that one of its 
ends is fixed in a cross section of a long cone-shaped 
channel and use the blob method to calculate the ef- 
fective driving force. For a conical channel, the diam- 
eter of a blob depends on its position along the chan- 
nel. It is equal to the channel diameter at each position; 
£(x) = Z?o + 2(a + x) tana (see Fig. [5]). Here, a is the 
distance of the fixed end of the polymer from a cross 
section of the channel with diameter Do. Following Ref. 
fl8j , we write the confinement free energy of the polymer 

as fcfr ~ Jo W)> which S ives 



' r ln(Do + 2(a + L)tana) - In (D + 2a tan a)] . (3) 



ksT 2 tan a 

I 

Here, L is the polymer size along the channel. The num- Also, the number density of the monomers inside each 
ber of monomers inside each blob is g(x) ~ (^f^J " [l9j - 
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blob is p(x) ~ fMjs , and the cross-section area of each 

blob scales as £(x) 2 . Hence, the number of monomers 
inside a region of thickness dx inside the channel (the 
region colored in gray in Fig. [2]) can be written as 
n(x)dx ~ f(fy < ^ E 03 ■ The size of the polymer along the 
channel can be obtained by the equality of the integral of 
n{x) over x and the total number of the monomers, N; 



FIG. 2: Schematic of a polymer inside a long cone-shaped 
channel. One end of the polymer is considered to be fixed in 
a cross section of the channel. The blob diameter as defined 
in the first method of the theory (in the case of large apex 
angles) is shown in the figure. 
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The force exerted to the polymer can be calculated from 
the derivative of the confinement free energy, Eq. [3J with 
respect to a, at a = 0. Here, one should note that N is 
the total number of the monomers on the polymer and 
its value is constant. Instead, the polymer size along 
the channel axis, L, depends on the parameters such as 
a (Eq. [3J. From Eq. [3J we can calculate the force /; 
~r~Tp ~ - (l + (tt) n ) n xot ; 1" tt- Tne derivative 

KgT V Koala— _Do+2-Ltana Do 

of L with respect to a is found from Eq. 2J (§j) o _ ~ 
( g +2 J Ltan» )-" 1 " L Accordingly, 
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1 
~D~a 
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Dq + 2Ltana 
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To calculate the force in the polymer translocation 
through a cone-shaped channel, L in Eq. [5] should be 
substituted by the channel length. The assumption of 
constant N and variable L in the calculation of the force 
is also valid for the case of polymer translocation. Be- 
cause, at each moment, the polymer segment inside the 
channel does not feel the finite length of the channel. 

The exerted force can also be calculated in another 
way, whose algebraic calculations are more complicated. 
This re-calculation is useful for checking the results. For 



this end, the blobs are defined as spheres tangent to 
the channel wall that cannot penetrate into each other. 
The same method for defining the blobs has been used 
in Ref. [2(J. The blobs are assumed as spheres tan- 
gent to the cone, so, the radius of a blob with its cen- 
ter at position x is = Do cos a + 2(x + a) sincv (see 
Fig. |3Ja)). Confinement free energy of the polymer is 
proportional to the number of these blobs. To count 
the blobs and calculate the confinement free energy, we 
use the geometrical relation between the positions of 
two consecutive blobs inside the channel; aij+i — x% — 
£(xi)+£(x i+1 ) ^ gee pjg^ [3Ja)). Using this recursive rela- 
tion, one can find the explicit relation between the po- 
sition of a blob, Xi, and its number along the channel, 
i; Xi = A l ~ 1 (xi + a + -2a cot a) — (a + ^& co t a) , where 



A = i+slna . The first and the last blobs are tangent 

1 — sin a ° 

to the beginning and the end of the channel, respec- 
tively. Hence, their positions along the channel, Xi and 



can be obtained from the relations 



and 



L, respectively, n is the total number of the 
blobs (see Fig. [3]). Using the equations for X{, X\ and 
x n , the number of the blobs, n, and the confinement free 
energy are obtained as 



J 



F log [2(L + a) sin a + Do cos a] — log [2a sin a + Do cos a] 

ksT log(l + sin a) — log(l — sin a) 

I 



(6) 



The polymer extension along the channel axis can be ob- 
tained from the constraint TV = J2i=i 9( x i)i where g(xi) 



is the number of monomers inside the zth blob. Substitut- 
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ing <7(xi) ~ (^f^-)" and^(a;i) ~ Dq cos a+2(xi+a) sin a, and using the equation for Xi, one obtains 



N 



[2(Z + a) sin a + D cos a] " — [2a sin a + D Q cos a] 
(1 + sina)£ — (1 — sina)^ 

I 



(7) 
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FIG. 3: (a) Definition of the blob in the second method of the 
theory (the case of large apex angles). The size of the blob at 
position x on the channel axis, and the relation between the 
positions of two consecutive blobs are shown in the figure, (b), 
(c) The first blob is tangent to the beginning of the channel, 
and the last blob is tangent to its end. 



Using the derivatives of Eqs. [B] and [71 we have 
/ 1 tan a 



k B T D 



o log (i 



1+sin at 
sin a 



l 



D 



Dq + 2L tan a 



(8) 

Note that Eq. [8] is quite similar to the result obtained in 
Eq. and differs only on the coefficient i^rpshri ■ 

When the diameter of the first blob is smaller than the 
length of the channel, £(xi) < L, we expect the results 
obtained from this method to be the same as those of the 
previous method. The value of £(xi) does not depend on 
the channel length, but it increases with the apex angle, 
rapidly. As is shown in Fig. |4j the two methods give 
the same dependence of the force on the channel length. 
Also, dependence of the force on the apex angle calcu- 
lated from the two methods is the same at small apex 
angles, and becomes different only at large angles. In 
the case of large apex angles, despite the case of small 
ones, the force increases monotonically with the apex an- 
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FIG. 4: (Color online) The force exerted to a polymer, which 
is confined inside an open cone-shaped channel, versus (a) the 
cone apex angle, with the channel length L — 15b, and (b) 
the channel length, with the apex angle a = 3°. In these fig- 
ures, Do = 2b. The force is calculated from the two methods 
described in the text. When the apex angle is small, the two 
methods have the same results. The force increases with the 
cone angle and length and becomes constant afterwards. 



gle (Fig. IUa)). The force increases with the channel 
length, L, and becomes constant, at larger values of L 
(Fig. IDJb)). Indeed, at large values of L, the channel 
diameter becomes larger than the radius of gyration of 
the polymer and it does not have any confining effect on 
the polymer. 

The length of the polymer segment inside the chan- 
nel, N, versus the channel length and the apex angle is 
shown in Figs. 03a) and[5jb). Although the values of N 
obtained from the two methods are different, they have 
the same dependence on the channel length and the apex 
angle. For comparison with the simulation results pre- 
sented in the next section, power-law functions are fitted 
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FIG. 5: (Color online) The length (number of monomers) of 
the polymer segment inside the channel versus (a) the channel 
apex angle (for L — 15b), and (b) the length of the channel 
(for a = 3°). To make the comparison easier, the polymer 
length obtained from the second method is multiplied by 3, 
in both panels (a) and (b). As can be seen, the two methods 
predict the same dependence of the polymer length on the 
channel length and the apex angle. Power-law fits to the 
curves are shown, for later comparison with the simulation 
results (see Sec. IIIip . 



where a is the monomer diameter and the MD length 
scale, and e is the Lennard- Jones energy scale. The 
monomers are connected by the harmonic potential, 
Uharm = ~^K(r~R) 2 , in which K is the spring constant, 
r is the distance between two consecutive monomers 
along the polymer and R is their equilibrium distance. 
The cone-shaped channel and the two walls are mod- 
eled by the Lennard- Jones potential between them and 
the monomers (see the schematic of the channel and the 
polymer in the inset of Fig. [6|). 

The time step of the simulations is r = O.Olro, where 

To = \J~^^- is the MD time scale, and m is the monomer 
mass. The simulations are performed in the NVT en- 
semble, using the Langevin thermostat, at the constant 
temperature T = le/ks- The Langevin equation mf — 
—VU(r) — £/,r + F ext + ff(t) is integrated, for describing 
the monomers motion, where £b is the friction coefficient, 
and F ext is the external force from the channel and the 
walls. ff(t) is the Gaussian white noise, which follows 
the fluctuation-dissipation theorem. The simulations are 
done with ESPResSo [H|. 

At the beginning of the simulation, the monomers are 
arranged on the axis of the channel, such that the lengths 
of the two segments of the polymer, which are outside 
the channel from the two sides, are the same. Then, we 
fix the monomer in the channel apex and let the other 
monomers to equilibrate. After that, we release the fixed 
monomer and set the simulation time equal to zero. The 
passage time is the time that the last monomer of the 
polymer lies in the channel apex. In our simulations, 
the diameter of the channel apex is Dq = 2.46 and the 
polymer has 100 monomers, unless otherwise stated. 



B. Simulation results 



to the results. In Fig. [Ha), two power- law functions 
with exponents 0.20 and 0.50 are fitted to the result of 
the theory, at small and large apex angles, respectively. 
Also, as is shown in Fig. [5jb), the curve obtained from 
the theory is followed well by a power-law function with 
exponent 1.15. 



III. SIMULATION 

A. The simulation method 

We use molecular dynamics (MD) simulations to check 
our theory. We use the coarse-grained bead-spring model 
for the polymer. The interaction between the monomers 
is the short-ranged Lennard- Jones repulsive potential 

ULj[r) ~ \ r > 2 1 /6 CT; 

(9) 



When the channel is cylindrical (a = 0), the poly- 
mer exits the channel with the same probabilities from 
the two sides of the channel. The polymer passage time 
in this case scales with the polymer length by the ex- 
ponent 2.2. This is consistent with the previous results 
for unbiased polymer translocation (22j. Instead, in the 
translocation from the cone-shaped channel, the poly- 
mer leaves the channel from the base in all our sim- 
ulations and the passage time scales with the polymer 
length with the exponent 1.4 (see Fig. [5]). This exponent 
is close to the exponent 1.6 predicted in the literature, 
for the driven translocation (of infinitely long polymers) 
[l5l [23j . It also completely matches the exponent for the 
driven translocation of shorter polymers [22| ■ This shows 
that the force exerted from the cone-shaped channel is 
determining and the polymer translocation through this 
channel is a driven process. 

The polymer passage time versus the channel apex an- 
gle is shown in Fig. [TJa), for L — 156. Upon increasing 
the channel angle, a, from °, the passage time decreases 
and shows a minimum at 2 . Then the passage time 
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FIG. 6: (Color online) The translocation time of a polymer 
through a cone-shaped channel, versus the polymer length. 
The channel parameters are L = 156 and a — 3 °. The expo- 
nent is close to the case of the driven polymer translocation 
[23| . The error-bars are resulted from averaging over 21 runs 
of the simulation. Inset: schematic of the cone-shaped chan- 
nel, the two walls and the polymer used in the simulation. 



increases up to its maximum value at 25 ° and then de- 
creases at larger angles. Note that at small angles, the 
channel base diameter is small and the simulation data 
should be described by the case of small apex angles of 
the theory, Eq. [2] The force obtained from this equation, 
first increases and then decreases with the apex angle, as 
shown in Fig. [1(a). Considering that the passage time 
is a decreasing function of the force, it would have the 
same behavior as the simulation result. However, at large 
apex angles, the channel base diameter is large and Eq. [5] 
should be used to describe the force exerted to the poly- 
mer. The force that is obtained from the case of large 
apex angles of the theory increases monotonically with 
the angle (Fig. Ufa)). This is in agreement with the 
passage time reduction at large apex angles. Combina- 
tion of the predictions of the two cases of small and large 
apex angles of the theory describes the simulation result, 
reasonably well. 

In the theory, in the case of the small apex angles, we 
assumed that the channel is a closed volume. To justify 
this assumption, the average density of the monomers in 
the range of distance 3cr from the channel base is calcu- 
lated in the simulations (see the inset of Fig. [7(a)). It 
can be seen that the monomers density close to the chan- 
nel base is higher, in the case of small apex angles. The 
sudden decrease in the value of the monomers density 
occurs where the base diameter becomes of the order of 
the monomer diameter; L tan(o) ~ b. 

The polymer passage time versus the channel length is 
shown in Fig. [7(b) , for the apex angle a — 3 ° . This time 
decreases with the channel length and then becomes con- 
stant. From the theory, the force in the case of small base 
diameters increases with the channel length, at small val- 
ues of L (Eq. [2] and Fig. H(b)). Also, in the case of large 
base diameters, the force increases with L monotonically 
(Eq. [5] and Fig. E(b)). These are in agreement with the 
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FIG. 7: (Color online) The simulation results for the passage 
time of the polymer versus (a) the channel apex angle (for 
L = 156), and (b) the channel length (for a — 3°). The 
passage time is a non-monotonic function of the apex angle, 
but decreases monotonically with the length. The inset of 
panel (a) shows the monomers density outside the channel 
base. It shows that the monomers crowd close to the channel. 
The error-bars are resulted from averaging over 15 runs of the 
simulation. For the zero apex angle, 42 runs of the simulation 
are done. 



result shown in Fig. [7(b). 

It is worth studying the number of monomers of the 
polymer segment inside the channel during its passage 
through the channel. The mean value of the monomers 
inside the channel as a function of the channel angle and 
length is shown in Figs. [5(a) and [5(b). The number of 
monomers versus the channel length can be described by 
a power-law function with exponent 1.14. However, the 
number of monomers against the tangent of the chan- 
nel apex angle can be fitted with two different exponents 
0.13 and 0.38, for the small and large apex angles, respec- 
tively. These figures are plotted for the same parameters 
as those of Fig. [5] and the fit exponents are nearly close. 



IV. SUMMARY AND DISCUSSION 

An important assumption in our theoretical calcula- 
tion of the entropic force exerted to the polymer is that 
the polymer is in equilibrium inside the channel. How- 
ever, it is known that the polymer passage through a 
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£f, the friction constant for each monomer. In this rela- 
tion, the friction force acting on the polymer is calculated 
from the Rouse model [23 |. Therefore, the time needed 
for the polymer to traverse the channel length would be 
tl ~ ~ LN J^ b [l5j . The relaxation time of the poly- 
mer segment inside the channel from the Rouse model 
is r eq ~ TbN 1+2u , where N is the length of the poly- 
mer segment inside the channel. Here, one should note 
that by using such scaling relations, that their numerical 
pre-factors are not known, judgment on our equilibrium 
assumption is not possible. For this reason, we perform 
simulations to check the polymer equilibrium inside the 
channel. In separate simulations, one end of a polymer 
is kept fixed at the channel apex and after equilibration 
of the polymer, the density profile of the monomers in- 
side the channel is measured . The density profile of 
the monomers inside the channel is also obtained during 
the passage of a long polymer (in the simulations of Sec. 
IIII[) . Comparison of the results of the two simulations 
shows the assumption that the polymer segment inside 
the channel is in equilibrium, is reasonable, in the range 
of the parameters used in our studies. This shows that 
in our simulations the polymer length inside the channel 
is small enough to equilibrate in a short time. 



FIG. 8: (Color online) The number of the monomers of the 
polymer segment inside the channel as a function of (a) the 
channel apex angle (for L = 156), and (b) the channel length 
(for a = 3 °). As is shown, the fit exponents are close to those 
of Fig. El 

channel is a non-equilibrium process, and it has been 
shown that the two ends of the polymer out of the chan- 
nel cannot equilibrate during the passage process [l5l[T6| . 
Here, we investigate the validity of the assumption that 
the polymer segment inside the channel is in equilibrium 
during the passage process. For this end, we compare the 
relaxation time of the this segment with the time needed 
for the polymer to traverse the channel length, L. If the 
former time is smaller, the assumption is reasonable. The 
velocity of the polymer passing the channel can be writ- 
ten as v ~ -Aft ) where / is the entropic force exerted 
to the polymer, N p the total length of the polymer, and 



In summary, the effective force of entropic origin act- 
ing on a polymer in the course of its passage through a 
cone-shaped channel was calculated and compared with 
the results of coarse-grained MD simulations. The force 
was obtained for the two cases of small and large apex 
angles of the channel. Combination of the results of the 
two cases showed that the effective force exerted to the 
polymer is a non-monotonic function of the channel apex 
angle. The force increases monotonically with the chan- 
nel length. The simulations showed the importance of 
the force exerted by the channel. It was shown that the 
simulation results for the polymer passage time through 
a cone-shaped channel can be described with the theory. 
Also, it was shown that the simulation and the theoretical 
results for the polymer length inside the channel during 
the passage of a long polymer are in good agreement. 
The simulation results also support the assumption that 
the polymer segment inside the channel is in equilibrium 
during the polymer passage. 
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